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Abstract 

Optimizing resonator shapes for maximizing the ratio 
of maximum to minimum gas pressure at an end of the 
resonator is investigated numerically. It is well known 
that the resonant frequencies and the nonlinear standing 
waveform in an acoustical resonator strongly depend on 
the resonator geometry. A quasi-Newton type scheme 
was used to find optimized axisymmetric resonator 
shapes achieving the maximum pressure compression 
ratio with an acceleration of constant amplitude. The 
acoustical field was solved using a one-dimensional 
model, and the resonance frequency shift and hysteresis 
effects were obtained through an automation scheme 
based on continuation method. Results are presented for 
optimizing three types of geometry: a cone, a horn-cone 
and a half cosine-shape. For each type, different 
optimized shapes were found when starting with 
different initial guesses. Further, the one-dimensional 
model was modified to study the effect of an 
axisymmetric central blockage on the nonlinear 
standing wave. 

Introduction 

The waveform of the standing wave in an oscillating 
closed cavity is strongly influenced by the geometry of 
the cavity (resonator). Lawrenson et al . 1 at MacroSonix 
Corp first exploited the dependence for developing 
high-amplitude acoustic pressures, referred to as 
resonant macrosonic synthesis (RMS). Peak acoustic 
pressures that exceed the ambient pressure three to four 
times and ratios of maximum to minimum pressures of 
27 were observed in RMS cavities. The size of the 
demonstrated overpressure reached the level that 
required by commercial applications such as acoustic 


pump or compressor. They considered these types of 
axisymmetric shape: cylinder, cone, hone-cone and 
bulb, and appeared to find that the hone -cone shape 
generates higher overpressure given an input power. A 
companion paper by Ilinskii et air developed a one- 
dimensional time domain model for studying the RMS 
numerically. Their results confirmed the nonlinear 
standing waveform and the related characteristics such 
as shape-induced resonance hardening and softening 
observed in the experiments by Lawrenson et al . 1 Chun 
and Kinr numerically investigated the cosine shapes in 
addition to cylindrical and conical shapes using high- 
order finite-difference approximations. They claimed 
that the half cosine-shape is more suitable to induce 
high compression ratio than other shapes. To serve 
commercial needs, such as in acoustic gas compressor 
and acoustic liquid pump, the common goal in these 
studies is to find an optimized shape for generating 
higher overpressure. Flowever, the optimization 
procedure and result have not yet been discussed. 

In this article, we introduce the numerical schemes and 
present the results in searching for the optimal shape 
parameters that give us the highest maximum-to- 
minimum pressure ratio in each of the following types 
of shapes: cone, horn-cone and half cosine-shape. In the 
second part, we study the effect of a cylindrical central 
blockage on the RMS. 

I. Governing Equations 

In this section, we describe a one -dimensional model 
for the acoustic wave field in the resonator that includes 
a central blockage. We develop our model following 
the work by Ilinksii et al. 2 for an axisymmetric 
resonator in absence of the central blockage. 
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Consider the acoustic field in an oscillating resonator 
driven by an external force. As sketched in Fig. 1, the 
resonator is axisymmetric with the inner radius given 
by r=r(x), 0<x<l, where x is the coordinate along the 
axis of symmetry. The central blockage is also 
axisymmetric with the same axis of symmetry x, whose 
shape is written as K = f 0 (x) . 

The density of the gas p, the velocity u and the pressure 
p satisfy the conservation of mass 


dp 

~dt 




~ r o)P U ) = °» 


( 1 ) 


and the conservation of momentum 


du du 1 dp . 
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dt dx p dx 
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a k e and 

a H = a ki(^ k ,x)J k = f k (v k ,^ k ,x). The 
detailed expressions for a k! and f k are the same as those 
in Ilinksii et air except that the central blockage is 
taken into account. 
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where o = , C 0 is the small signal 

Po 

propagation speed, D' m is given by 


where a(t) is the acceleration of the resonator enforced 
by the external force; ^ and r| are coefficients of 
viscosity. The state equation is specified by that of an 
ideal gas 


D' m = ~(r - 1 )xd m - im(y - X)dp m 

2 —— — [v 2 ] 

2{r 2 -r 2 f V im ’ 


P = Po{p/PoY . 


(3) and 


where p 0 and p 0 are the ambient pressure and density 

respectively. The no-penetration boundary conditions at 
the two ends require that the velocity vanish at x=0 and 
/. 

Following Ilinksii et al expressing the variables in 
finite Fourier series in time, the equations (1) and (2) 
can be reduced to a system of Ordinary Differential 
Equations (ODEs) for the Fourier coefficients of the 
velocity potential cp, defined as U = V <p : 


d(p k 

dx 

Tat, 

„ / 



dv, 

dx 


f k ’k = 
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where 


N 

[v 2 ]* = Zv/V,. 
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The no-slip boundary conditions at the two ends are 
translated to the Fourier space: V k = 0 at x=0,l. 


After the Fourier coefficients { (p k } and { V k } are 
obtained by solving the boundary-value problem Eq.(4), 
the velocity potential 9 and the modified velocity v are 
computed from the inverse FFT. The density p is given 
by the momentum equation 
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and the pressure p can be obtained from the state 
equation p / p 0 = (/? / p 0 ) r . 


of hysteresis effects). The method for obtaining Re for 
fixed resonator shape is explained later in the section. 


II. Numerical Methods 

In this section, we describe the procedures for finding 
the optimal shape parameters so that the desired 
standing waveform can be obtained. 

To simplify our discussion, we introduce the following 
dimensionless variables 


X=~, T = cot, R=~, 

l l 




A =fj’W 

lco n 


where co is the frequency of the periodic force acted on 

the resonator, co 0 = — ~ — is the fundamental frequency 

of a cylindrical resonator of length / . 

In this work, we assume the acceleration of the 
resonator is harmonic, A(T) = Acos(T) . The 
acoustic wave field in the resonator is determined by 
the acceleration amplitude A , the ratio of specific 
heats y, the attenuation coefficient G, defined as 

G = + _ anc j q le resona ]; 0r oscillating 

c oPo 

frequency Q, defined as Q = — . 

©o 

Now, we briefly describe the method for finding the 
optimal resonator shape. 

In the optimization process, we fix the values of the 
acceleration amplitude A = 5 x 1(T 4 , the gas specific 
heat ratio y=1.2 and the viscosity -related parameter 
G=0.01. Suppose the resonator shape R(X) is 
determined by a number of shape parameters, S 0 , S l5 
. . . , S n For example, a cone is given by 
R{X) = S 0 + S^X . The compression ratio of a 

resonator R c , defined at the ratio of maximum pressure 
to the minimum pressure at the narrow end of the 
resonator 


KAoA at.V 0. 


( 6 ) 


is a function of the shape parameters, the dimensionless 
frequency Q and the history of Q (due to the existence 


We first outline the optimization schemes. A quasi- 
Newton method, BFGS (Broyden-Fletcher-Goldfarb- 
Shanno 4 ), is used for maximizing the multi-variable 
nonlinear function R c (S 0 ,S l S ) . Since the 

evaluation of the objective function Rc itself involves 
solving a nonlinear system of ODEs Eq. (4) for many 
times, the gradient information of Rc required for the 
BFGS method is not available analytically and is 
derived by partial derivatives using a numerical 
differentiation method via finite differences. This 
entails perturbing each design variables, S h in turn and 
calculating the rate of change in the objective 
function R c {S 0 ,S l S n ) . For two shape parameters, 

the optimization takes 4 to 48 hours of CPU time 
(depending on the type of the resonator shape and 
initial guess for the shape parameters) on a 1 .3 GHz 
Athlon T-Bird PC with the Lahey-Fujitsu FORTRAN 
compiler. 

For a given shape of the resonator and the central 
blockage, the boundary value problem Eq. (4) is solved 
numerically by a Multiple Shooting Method. Because 
of the hysteresis effects, the solution is not unique near 
the resonant frequency and the Multiple Shooting 
Method will not converge unless a good initial guess of 
the solution is provided. To circumvent the difficulty, a 
continuation method is implemented: the system of 
ODEs is solved starting from a frequency Q that is far 
away from the resonant frequency and the solution is 
used as an initial guess for solving the ODEs for 
increased or decreased Q; the steps are repeated until 
all branches of the solution for all values of Q near the 
resonance is completed. The maximum ratio R c among 
different values of Q is chosen as the compression ratio 
for the resonator. 


III. Results 

A. Characteristics of the standing waves 

Before we present our results on finding the optimal 
resonator's shape, we illustrate some of the important 
properties of the pressure wave in a non-cylindrical 
resonator, whose shape is described by that of a 
homcone: 

js 0 coshOVO,0 < X < S 2 
\a + j3X,S 2 <X<\ 
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where 

S 0 =0.028333,5; =5.7264,5', =0.25,a = 5' 0 cosh(5' 1 5' 2 ), 
and p = S 0 sinh(<S' 1 <S’ 2 ). The shape parameters for 

the homcone are obtained from those used in tire 
experiments by Lawrenson et al. 1 

In Fig. 2, we show the pressure waveforms at the ends 
of cylindrical resonator and the horncone resonator at 
their corresponding resonance frequencies for the 
effective viscosity G=0.01, the gas specific heat ratio 

y=1.2, the acceleration ^4=10 3 . Throughout this 
paper, we will keep the effective viscosity G and the 
gas specific constant at these values, which are close to 
the values used in the experiments 1 and our numerical 
simulations have suggested that the pressure waveform 
is not a strong function of these parameters. At the 
same level of the acceleration, the waveform at the 
narrow end of the horncone shows large variation in 
which the ratio of the maximum and minimum 
pressure, R c , exceeds the value of 30.5. On the other 

hand, the waveform for the cylinder shows the 
formation of shocks at the resonance and the 
compression ratio R c is below 1.32. The difference in 

the waveforms shows the strong dependence of the 
acoustic field on the resonator shapes. The graph also 
shows that the variation in pressure at the wide end of 
the homcone is much milder than that of the narrow 
end, oscillating within 23% above or below the value of 
the reference pressure p 0 . 

In Fig. 3, the ratio of the amplitudes of the second and 
the first harmonics of the pressure, /A / p x , are plotted 
against the reduced amplitude of the first harmonic. 

Pi / P 0 , calculated at one end of the cylinder and at 

the narrow end of the homcone with increasing 
frequency. The figure shows that the second harmonic 
reaches its maximum at small amplitude of the first 
harmonic for the consonant (cylinder) resonator and the 
second harmonic increases slowly as the first harmonic 
increases for the dissonant (homcone) resonator. In 
other words, the energy is pumped into the first 
harmonic more for the horncone than that for the 
cylinder, creating the large-amplitude pressure wave in 
the homcone resonator. In contrast, the energy is 
efficiently absorbed by all harmonics of the pressure 
wave in a cylindrical resonator and the amplitudes of 
the harmonics reach their maximums at small amplitude 
of the pressure. 

In Fig. 4, the amplitude of the first harmonic is plotted 
against the frequency of the oscillating resonator for 
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two different levels of acceleration: A — 2 X 1(T 4 and 

A = 1(T ! . For the smaller acceleration, the pressure is 
uniquely determined at each frequency; for the larger 
drive amplitude, the pressure takes one of the multiple 
values near the resonance, depending on the direction 
of the change in frequency. The existence of hysteresis 
and hardening resonance in homcone resonators 
requires that the largest pressure compression ratio be 
obtained through an upward frequency sweep. Our 
numerical simulations show that the frequency 
increment size near the resonance must be small for the 
convergence of the Multiple Shooting method. 

B. Optimal conical shapes 

In following subsections, we present results on 
optimizing resonator shapes to achieve maximum 
compression ratio R c at one end of the resonators. For 
simplicity and reducing computing cost, we assume the 
central blockage is absent, i.e., R 0 =0, and keep ratio of 
specific heats y=l .2, the attenuation coefficient G=0.01 
and the acceleration A = 5 x 1(T J constant. 

A conical resonator can be written as 

R(X) = S 0 + S ] X, for 0<X<1. (8) 

Using the optimization procedure described in Sec. II 
and starting with the conical shape given in Lawrenson 
et al ., 1 we found that the compression ratio R c of a 
conical resonator reaches a maximum when the cone 
has the shape constants: (S^O. 032945 and Si=0. 26800. 
The achieved compression ratio is R c =5.048 for the 
fixed values of the parameters given earlier. Recall that, 
due to hysteresis, in order to find R c> one has to trace 
out the entire branch by incrementing the frequency Q. 
The compression ratio reaches the value when the 
frequency is increased to 1.3135. The conical resonator 
studied in Lawrenson et al . 1 and Ilinksii et al . 2 has 
almost identical shape Stf=0. 032941 and .S^O. 26800. 

The adopted optimization scheme BFGS is designed for 
finding a local extreme of a multivariable function. Our 
numerical simulations indicate that the compression 
ratio, as a function of the shape parameters, have many 
local extrema. Starting with a different initial shape, we 
obtained another locally optimized conical shape with 
S 0 =0. 010413 and Si=0. 19572. When the frequency is 
raised to the value of 1.3862 by small increments, the 
compression ratio reaches the value of 5.343, which is 
about 6% higher than that for the first optimal conical 
resonator. In the current optimization, the slope of the 
cone is restricted to be less than 0.268 and we require 
that the radius at the small end be larger than 0.01. The 



two optimal conical resonators are shown in Fig. 5. The 
second cone (Fig. 5(b)) has the smaller narrow end and 
smaller slope than the first one (Fig. 5(a)). 



Horn-cone geometry is described by 


R(X) 


J R h cosh(D h X) for 0 < X < X c 
\a+/3(X-X c ) for X C <X<\ 


frequency Q=1.5683, which is about 1 1% better than 
that of Chun and Kim 3 at the same level of acceleration. 
Under the specified conditions, we found that horn- 
cone shape is better than the cosine shape in generating 
higher compression ratio at the narrow end, which is 
opposite to the findings of Chun and Kim. 3 Shown 
in Fig. 7(c), the pressure waveform of the optimal 
14-cosine design is smoother than that of Chun and 
Kim. 3 


E. The effect of the central blockag 


where 

a = R h cosh (D h X c ) and J3 = R h D h sinh (D h X c ) . 

Rh, X c and Dh are positive constants. We consider that 
the connection point of horn and cone in the homcone 
is fixed atX,, =0.25. We optimize the remaining two 
shape parameters S 0 = R h and S x = D h . We found 

that the dimensions of the horn-cone given in 
Lawrenson et al.\ Rh =0.028333, Dh =5.7264, are 
nearly locally optimized if we start the optimization 
with these values. For these parameters, the maximum 
compression ratio 1 3 . 1 0 is achieved when the frequency 
is increased to 1.4671. In searching the optimal 
horncone design, the parameters, such as acceleration 
and specific heat ratio, are fixed at the values given in 
the beginning of Sec. Ill B. Starting with a different 
initial geometry of the horncone, we can obtain a higher 
value of the compression ratio R c =\l A2 when R h 
=0.024488, Dh =5. 1434. For the compression ratio, the 
frequency was increased to 1.41 10 gradually with small 
increments. Comparing with the compression ratio for 
the homcone in Lawrenson et al . 1 (shown in Fig. 6(a)), 
the second optimal homcone (shown in Fig. 6(b)) 
improved the ratio by more than 30% at the same value 
of acceleration. Shown in Fig. 6(c), the pressure 
waveform at the narrow end for new homcone design is 
a little more complicated than that for the original 
horncone given Lawrenson et al . 1 



We define a 14-cosine shape as 

R(X) = S 0 + Sfl-cos(nX), for 0<X<1. 

The parameters of the resulting optimal cosine shape 
are S 0 =0.0 16307, £,=0.093041, shown in Fig. 7(b), 
compared with £<,=0.025, £,=0.095, the !4-cosine 
dimensions reported in Chun and Kim. 3 For the optimal 
shape, we obtain the compression ratio Rc =10.67 at the 


In this section, we discuss the influence of the central 
blockage on the pressure waveforms in a resonator of 
fixed shape. In this work, we limit ourselves to 
cylindrical central blockage, though our equations 
derived in Sec. I apply to arbitrary axisymmetric shape 
of central blockage. 

Figure 8 shows the pressure waves at the narrow end of 
a conical resonator described by Eq. (8) with 
£<f=0. 032945 and £,=0. 26800 for different central 
blockage radius sizes: R 0 = 0, 0.01 and 0.02. These 
pressure profiles were obtained at the fixed frequency 

Q = 1.3132 and for A = 5 x ] 0 4 , G = 0.01, y = 1.2. 
The central blockage in this case was a simple cylinder. 
It appears that the presence of a central blockage 
reduces the amplitude of the pressure wave. However, 
Fig. 9 shows the resonant frequency of the conical 
resonator with a central blockage is shifted. Keeping 
other parameters constant, the pressure waveforms for 
the conical resonator with the cylindrical blockage of 
radius R 0 = 0.02, at different oscillating frequencies Q = 
1.2812, 1.3068 and 1.3291, are shown inFig. 9. 
Comparing to the pressure wave shown in Fig. 8 for a 
resonator without a central blockage, the amplitude of 
the pressure is nearly the same when the central 
blockage of radius 0.02 is included. The results indicate 
that a cylindrical central blockage in a conical resonator 
would affect the resonant frequency but would have 
negligible effect on the amplitude and overall shape of 
the pressure wave at the resonant frequency of the 
resonator/blockage system. 

IV. Conclusion 

We have performed local optimization schemes for 
finding the resonator shapes that maximize the pressure 
compression ratio at one end of the resonators. The 
dimensions for several types of optimal resonators are 
reported: cone, homcone and 14-cosine shapes. For each 
type, we found there are many different designs that 
achieve local extrema. Trying with different initial 
guesses for the optimal design, we show that one can 
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get as much as 30% improvement on the compression 
ratio with a fixed level of acceleration. For the shapes 
we considered, it appears that the horncone shape 
generates the highest compression ratio. Strategies for 
searching globally optimal shapes are imder 
investigation. 

The effect of including a cylindrical central blockage in 
a resonator is studied in this work as well. For a conical 
resonator, it seems that the resonance frequency is 
increased when a central blockage is present and the 
compression ratio at resonance is not significantly 
affected. 
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Figure 1. Sketch of an axisymmetric resonator and a central blockage. 



Figure 2. The reduced pressure p / p 0 at the ends of a horncone and a cylinder at their corresponding 
resonance frequencies for the same acceleration A = 10 
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Figure 3. The ratio between the amplitudes of the second and the first harmonics of the pressure wave, p 2 / p l , 
is plotted as a function of the dimensionless amplitude of the first harmonic reduced pressure p x l p 0 . With 
increasing frequency, these amplitudes are calculated for the narrow end of the horncone and at one end 

of the cylinder. 



1.35 1.4 1.45 1;,5 


Q 

Figure 4. The effect of the frequency and the amplitude of acceleration on the dimensionless amplitude of the 
first harmonic reduced pressure p x / p 0 . The pressure is calculated for the narrow end of the horncone for the 

acceleration A = 2x10 4 (crosses) and A = 1 0 3 (circles and squares). 
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DC 


DC 


Figure 5. Optimized conical shapes: the compression ratios corresponding to the shapes in (a) and (b) reach the 
values of 5.048 and 5.343 respectively. For shape parameters, see the text. 






Figure 6. (a) The homcone shape in Lawrenson et al . 1 (b) Optimized homcone shape. For shape parameters, 
see the text, (c) The pressure waveform at the narrow end for the homcone in Lawrenson et al . 1 (the solid line) 
and that for the second optimized homcone (the dashed line) are shown. 
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(a) (b) 





Figure 7. (a) The i/ 2 -cosine shape as in Chun and Kim 3 , (b) The optimized '/ 2 -cosine shape, (c) The pressure 
waveform at the narrow end for the '/ 2 -cosine resonator in Chun and Kim 3 (the solid line) and that for the 
optimized '/ 2 -cosine shape (the dashed line) are shown. 



Figure 8. Effect of the center blockage. The pressure waveforms at the narrow end of a cone resonator are 
plotted for the different radii of the central blockage: R 0 = 0 (solid line), 7? o =0.01 (dashed line) andi? 0 =0.02 
(dotted line). The parameters corresponding to the plot are given by Q = 1.3132, A = 5 x 1(T 4 , G = 0.01, y = 1.2. 
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Figure 9 . Pressure waveform at different frequencies: Q = 1.2812 (the dotted line), Q = 1 .3068 (the dashed line) 
and Q = 1 .3291 (the solid line), for a conical resonator with the cylindrical central blockage of radius R 0 = 0.02. 
The rest of the parameters corresponding to the plot are given by xl = 5 X 1 0 , G = 0.01 , y = 1 . 2 . 
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